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We study the first derivative of tfie staggered magnetization squared dm^{9)^/d6 and the second 
derivative d^eo{9)/d9^ of the ground state energy per site. The parameter 9 controls the anisotropy 
between horizontal and vertical couplings in a two-dimensional (2D) spin-1/2 antiferromagnetic 
Heisenberg model. It is shown, that both derivatives diverge at 9 = 1, where the anisotropic 2D 
model reduces to the ID model. 



I. INTRODUCTION 

Antiferromagnetism looks different in the ID and 2D 
Heisenberg model with nearest-neighbor coupling of spin- 
1/2 matrices. The zero-temperature staggered magneti- 
zation is known t&,be zero in the ID and to be non- 
zero in the 2D model.cTEl Introducing different strengths 
for the horizontal and vertical coupHngs = J{1 + 0)/2 
and Jt, = J(l — 9)/2 the staggered magnetization m'^{d) 
has been studied as a function of the anisotropy pa- 
rameter 9. It has been found in Ref. ^ m^9)^ for 
< 6 < do, where 9a is at least Oq < 0.98E It has 
been suggested, that 771^(6*) ^ in the whole interval 
< 6* < 1. 

In this paper we study the derivative dm\6) /dO and 
the second derivative d'^eo{9)/d9^ of the ground state en- 
ergy. The transition from one to two dimensions is per- 
formed with an interpolating Hamiltonian H(0), defined 
on system sizes N = k'^ zL 1, k = 3,5,... with helical 
boundary conditions. The 0-evolution of energy eigen- 
values and eigenvectors as well as expectation values of 
hermitian operators is discussed in Sec. II. These gen- 
eral results are applied and numerically evaluated for 
dm^{9f/d9 and d^eo{9)/d9^ in Sec. Ill and IV, respec- 
tively. Our main result states that both derivatives di- 
verge logarithmically with the system size iV at 6* = 1 . 

The origin of these divergences can be traced back to 
the low energy excitations in the sector with total spin 
of the ID antiferromagnetic Heisenberg model. This sec- 
tor is investigated in Sec. V. In Sec. VI we discuss con- 
sequences for the staggered magnetization in the ID an- 
tiferromagnetic Heisenberg model with next to nearest- 
neighbor couplings. 



and N = k'^ zL 1 sites. The interpolating Hamiltonian 
between one and two dimensions is defined by: 

H(0) ^ ii^Hi + i^QfcHiQfct. (2.2) 

Qk represents a permutation operator for the sites on the 
ring 

QkS:Ml = S^^_,),+„ a = x,y,z, (2.3) 

which transforms the nearest-neighbor couplings in 
Hamiltonian (2.1) into couplings between spins separated 
by fc = VNtI- 



N 



QfcHiQ;, 



t _ 



(2.4) 



n=l 



H(0) defin es , a 2D Hamiltonian with helical boundary 
conditions, MO as can be seen in Fig. 1 of Ref. |l^. The 
horizontal and vertical couplings have strengths {l + 9)/2 
and (1 — 6')/2, respectively. The parameter 9 controls the 
interpolation from the ID Hamiltonian H(l) = Hi to the 
2D Hamiltonian with equal strength for horizontal and 
vertical couplings 



H4 



H(0), 



where 



H± = i(Hi±Q,HiQi 



(2.5) 



(2.6) 



The permutation operator is essentially a rotation, 
which transforms horizontal nearest-neighbor couplings 
into vertical ones. One verifies the equations: 



II. THE INTERPOLATING HAMILTONIAN 

Let us start with the ID Heisenberg Hamiltonian with 
periodic boundary conditions: 



N 

Hi = S„ • Sn+i, 



Sjv+1 — Si, 



(2.1) 



ca 
ca 



P = N +1 

P = N-1 



(2.7a) 
(2.7b) 



Therefore, Hi commutes with Q^, but it docs not com- 
mute with Qfe: 



[QLHi]=0, [Qfc,Hi]^0. 



(2.8) 



1 



The oper ator transforms the interpolating Hamilto- 
nian (2_^) as follows: 



Note, that is a unitary operator 



QfeQ 



(2.9) 



(2.10) 



and for this reason H{0) and H(— 0) are unitary equiva- 
lent and possess an identical spectrum of eigenvalues: 



H(0)|*„(0)) =S„(0)|*„(0)), 
Qfc|*«(e)) = |*„(-^)). 



(2.11a) 
(2.11b) 
(2.11c) 



Let us investigate the 6'-evolution of the eigenvalues 
En{0 ) and the eigenvectors |5'„(6')). Differentiation of 
Eq. (2.11a) with respect to 6 yields: 



d9 



Smn - [ErniO) - (0)] (^^) | (0)) 

+ (VI/„(0)|H_|\I'„(0)). (2.12) 



For n = m we find, that the (^-evolution of the energy 
eigenvalues 



dEn{9) 



1 

W 



[En{9) - {^,,{0)\Ii{-0)\^J9))] , (2.13) 



is governed b y th e expectation values of the Q^-rotated 
Hamiltonian ( |2.9| ). The 0-evolution of the eigenvectors 
|*„(6')) reads 



de 



Em{0) — En{0) 



which follows from (2.12) for n ^ m. This yields for the 
^-dependence of matrix elements of hermitian operators 

o{e): 

^^{^„mo{e)\^n{e)) = (vi/,„(0)|o'(0)|vi/„(0)) 



Ek{e)-E„,{e) 



E^^^1NM^(*.W|0(^)|*.(^)). (2.15) 



Eii9)-E„i9) 



In particular we get for the ^-evolution of the expectation 
values ('Fto (6> ) |H(- 6>) | (0)) , which enter on the right- 
hand side of (2.13): 



-(*„,(0)|H(-0)|vl/,„(0)) 



\{^„{0)\ll{^0)\^.U9))\' 
En{0)-E^{0) 



(2.16) 



It shoul d be noted that the second term on the right-hand 
side of (|2T|) can be expressed in terms of the dynamical 
structure factor: 

s^[u,-u{-0)]^Yl I(*«(^)|h(-0)|*™(0))|' 

x6{Lu-[E,n{,0)-E,,{9)]), (2.17) 
associated with the transition operator Yi{—0). 



III. THE FIRST DERIVATIVE OF THE 
STAGGERED MAGNETIZATION 



Equation ( 2.15 ) enables us to compute the slope of the 
staggered magnetization squared: 

mt(0)'^l(M.o(0)|St(vr)S(^)|vI/„(0)), (3.1) 

where 



(3.2) 



TO^(0)^ is known to be zero for 0= \ (D=l) and nonzero 
for 6* = (D=2) [ TOt(O) « 0.30 see Refs. It is 

expected that m^{0) > for < < 1. The derivative 
of the staggered magnetization is of special interest near 
6* = 1. Using ( |2T5| ) and (|3^) it follows: 



dwH0f_ 
dB~ 



M„,{0,n), 



LUm{0) 



(3.3) 



where we have introduced the notation: 



M^{0,p) = (vI/„(0)|S(-p) • 8(^)1*0(0)). (3.4) 

For the computation of the right-hand side of ( |3.3D we 
need the excitation energies LOrn{0) = Em{0) — Eq[0)\ the 
transition amplitudes (\E'm(6l)|H_|\l>o(6')), and Mra{0,ir) 
in the total spin S't = sector of the interpolating Hamil- 
tonian H(6'). Having determined the ground state vector 
\^q{0)) by means of the Lanczos algorithm, the excita- 
tion energies and transition amplitudes can be computed 
via the recursion method described in Appendix A. 

Going back to the definition ( |2.6| ) of H_ and express- 
ing the spi n operators S/ in terms of their Fourier trans- 
forms ( |3.2[ ), we arrive at the following representation of 
the transition matrix elements: 



(*„,(0)|H_|vI/o(0))=^ 



''P^)M„,{0,p). (3.5) 



It is important to note that the matrix elements Mm{0,p) 
are constrained by energy conservation for m ^ 0: 
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Q={<f^{9)\U{e)\^om, 



-f'^^ Mm{e,p). (3.6) 



The slope of the staggered magnetization ( p.3\ ) can now 
be written as: 



d0~ 



I.7T-27V/N , 

2/ -5[cos(pfc) - cos]5]E(6l,p,7r,iV), 

JO 



TT 



where we have defined: 



M„,i0,Pl)Mm{O,P2) 
UJm{0) 



(3.7) 



(3.8) 



S(6>j ^p, TT, A^) has to obey the sum rule resulting from 

-E(0,7r,^,iV) - / J^E{e,p,n,N) 

x[(l + 6i)cosp+ (1 - 6i)cos(pfc)]. (3.9) 

We first discuss the situation for 9—1 (D=l). As can 
be seen from Fig. 0, E(l,p, 7r,iV) is negative for 7t/2 < 
p < TT, whereas S(l, tt, tt, iV) is positive by definition. 
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FIG. 1. E(l,p, 7r,7V) with the data of finite systems 
= 8, 10, . . . 28, (o). The solid dots (•) show the forecast 
i nferr ed from (3.14). The solid line represents the prediction 



(3.12), whereas the dashed li ne in the inset represe nts t he 
corresponding least square fit (3.11) to E(l, tt, tt, Af) (3.10). 



In the large N limit 
1 



-I](l,7r,^,iV)'^-=5°Aln — 



(3.10) 



diverges logarithmically with N, as is shown in the inset 
of Fig. |l]. From a least square fit we estimate: 



A = 0.0542(5), iVo = 1.089(5). 



(3.11) 



Owing to the sum rule (p.9|), the divergence (3.10) is 
related to the singularity 



E(l,p,7r,oo; 



p — >7r 



-A 



1 — p/tt 



(3.12) 



which we compare in Fig. 1 with finite system results for 
S(l,p, TT, A'^). Finite-size effects a re sm all for momentum 
values away from the sing ularity (|3T|). In the combined 
limit 

p-^TT, A^^oo, {l~p/Tr)N fixed, (3.13) 

we try to describe the finite-size effects with a finite-size 
scaling ansatz: 



E(l,p,^,A^) = S(l,p,7r,oo)G(z), 



(3.14) 



which we already used in Ref. |l2yi3| to analyse sin- 
gularities of static structure factors. The ansatz al- 
lows to predict I](l,p, tt, 2A^), for p = Tr(l — 4/A^) and 
2A^ ~ 32, 36, . . . , 56. This prediction is shown in Fig. 1 
by the dots (•). These points yield a smooth extrapola- 
tion of the numerical results for z = 4, A^ = 16, 18, . . . j_28 



which are already close to the predicted behavior ( 3.12 ), 
represented by the solid curve in Fig. |l|. The importance 
of the singularity (3.12) is obvious, it yields the lead- 
ing contribution to the derivative of m^O) in the limit 
A'' ^ oo: 



dTOt(6l)' 



de 



N- 



N- 



2A / — [cos(pfc) — cosp], 

J2n/N P 



I2tt/N 

-A In N. 



(3.15) 



We would like to stress, that the singularity (3.15) be- 
comes visible only by following the arguments presented 
above. In p articular, the sum rule (3.9) and the diver- 
gence ( 3.10| ) of the susceptibility I](l, tt, tt, A^) in the ID 



case are crucial for the derivation of ( 3.12| ). 

One might wonder whether the singularity of 
dm\e) /de at 61 = 1 can be seen as well in the 6- 
dependence. For this purpose we have evaluated ( |3.3| ) 
on systems of size A^ = fc^=Fl, fc = 3,5. The results for 
dm\e) /de as a function of 6* are shown in Fig. ||. 

For all systems the first derivative dm'' {e)^ / d6 is neg- 
ative and decreases in the range < e < em{N), passes 
a minimum at em{N) and then increases again, by pass- 
ing a turning point Ot^N) near = 1. How to reconcile 
this behavior with our bias for the thermodynamic limit 
of dm^e) /d91 We expect a decrease for < 6* < 1, a 
singularity 



dm''{ef 

de 



-oo, 



(3.16) 



if we approach the ID limit from the left and 
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FIG. 2. The der ivative drn' (9) /dS for finite systems, eval- 
uated via Eq. (3.3). 







(3.17) 



for 9 > 1. For 6 > 1 the horizontal and vertical couplings 
in (2.2) arc antiferromagnetic and ferromagnetic, respec- 
tively. In this situation, the staggered magnetization is 

supposed to vani sh. 

If the scenario ( |3.16 ) and ( |3.17| ) is correct, there should 
be a sharp discontinuity in dm'^ (9)'^ / d9 when we pass 
the transition point 9 — 1. For the finite systems - 
presented in Fig. ^ - this sharp discontinuity is washed 
out. The transition region starts already at the mini- 
mum 9m{N) < 1 and extends far beyond the transition 
point 9 — 1. Therefore, we expect that for increasing 
system size the minimum position 9m{N) converges to 
the transition point 9 — 1: 



whereas the minimum value diverges: 



dm) [9) 



de 



(3.18) 



(3.19) 



Comparing in Fig. || the results for iV = 8, 10, 24, 26, we 
see already that the slope of dm){9)'^ /d9 becomes steeper 
for 9 > 9ra{N). This feature could be interprete d as a sig- 
natur e for the emergence of the discontinuity (3.16) and 
(3.17) in the thermodynamic limit. However, it certainly 
requires larger systems to observe this process. 



IV. 



THE SECOND DERIVATIVE OF THE 
GROUND STATE ENERGY 



According to ( 2.13| ) the first derivative of the ground 
state energy per site eo{9) at 6* = 1, 



rfeo(g) 
d9 



:eo(l), 



(4.1) 



is given by the ground state energy eo(l) = —(In 2— 1/4) 
itself.Ej Here we hasf&jused the fact that the ground state 
expectation valuetJila 



1 

n' 



*o|H(-l)|*o) - (*o|SoSfc|*o) 



vanishes in the limit k = viV=pT — > oo. Combining 
( pl3| ) and ( pl6| ), wc find for the second derivative: 



d^eoi9) 
d9^ 



2 



E 



|(vI/,„(g)|H_|^o(g))| 

LOrniO) 



(4.2) 



The Fourier decomposition ( |3.5| ) of the transition matrix 
elements (5'„i(^?)|H_ |vI'o((?)) leads to the representation: 



d^eo{9) 
d9^ 



4 / —[coap-cos{pk)fE{9,p,p,N) 












TT TT 



X [cosp — cos(pfc)][cosp' — cos{p' k)]Yi{9 , p, p' , N). 



(4.3) 



Y,{9,pi,p2, N) is defined in (3.8) and numerical results 
for 9 — 1 and TV = 24 are given in Table |. Note, that 
S(6' = l,pi,p2, N) < for 7r/2 ^ P2 < Pi < tt, whereas 
J:{9 = l,p,p,N) > by definition. 



TABLE I. T,{e,pi,p2,N) [Eq. (|3^] for iV = 24 and 
Pi,P2=2-kI/N, Z = 4,5,...,12. 



6 



8 



9 



10 



11 



12 



0.005 

-0.001 
-0.001 
-0.002 
-0.006 
0.000 
0.002 
0.003 
0.002 



-0.001 
0.007 

-0.002 
-0.006 
-0.001 
0.001 
0.001 
0.001 
-0.001 



-0.001 
-0.002 
0.006 

-0.001 
0.000 
0.001 
0.001 
0.000 

-0.005 



-0.002 
-0.006 
-0.001 
0.018 
-0.000 
-0.000 
-0.000 
-0.002 
-0.011 



-0.006 
-0.001 
0.000 
-0.000 
0.030 
-0.002 
-0.003 
-0.005 
-0.022 



0.000 
0.001 
0.001 
-0.000 
-0.002 
0.050 
-0.007 
-0.012 
-0.045 



0.002 
0.001 
0.001 
-0.001 
-0.003 
-0.007 
0.096 
-0.029 
-0.106 



0.003 
0.001 
0.000 
-0.002 
-0.005 
-0.012 
-0.029 
0.231 
-0.371 



0.002 
-0.001 
-0.005 
-0.011 
-0.022 
-0.045 
-0.106 
-0.371 
1.051 



T,{9,p,p' , N) and T,{9,p,7T, N) are related via the en- 
ergy conservation relation (^.6[): 



[(1 + 61) cosp +(1-9) cos{pk)]E{9,p,p, N) 

-^{9,p,7:,N) = A{9,p,N), (4.4) 



where 



tt-2tt/N 



Ai9,p,N)^-— / dp'+ I dp'] 

\J0 Jj}+2tt/N J 

X [(1 + 61) cosp' + (1 - 61) cos(p'A:)]E(6l,p,p', N). 

(4.5) 

The sum rule ( |4.4| ) for = 1 yields important con- 
straints for p ^ TT, where S(l,p, tt, oo) develops the pole 



4 



( |3.12 ) depicted in Fig. This pole has to be compen- 
sated by a corresponding singularity in: 



S(l,p,p,oo) 



B 



1 — p/tt ' 



(4.6) 



which is indeed observable in the numerical data shown 
m Fig. I The finite -size effects can be treated again with 
a finite-size scaling ansatz of the type (3.14). The open 



symbols (o) represent the predictio n of this ansatz. The 
solid curve is a pole fit of the type (4.6) with: 



B = 0.053(10). 



(4.7) 
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FIG. 3. iV) for finite systems N = 8, 10, . .^8 (•) 

and the predicted pole (4.6) in conjunction with (4.7) for 
E(l,p,p, cx)) (solid line). 

Th e po le for p — > tt appears as well on the left-hand 
side (14) of the sum rule: 



A{9 = l,p,N) 



A-2B 
1 — p/tt 



(4.8) 



The right-hand side of ( [4.5D tells us, how '£{9^1, p, p', N) 
behaves for p, p' — s- tt in order t o ge nerate the pole term 
(18) on the right-hand side of (4.4): 



N\nNJ:{l,p,p\N) 



(i-pA)(i-p'A)- 



Insertion of the leading singularities 
into the right-hand side of (O) yields 



and (4.9) 



de 



*=1 



A^A + 2B)\nN. 



(4.10) 



In Fig. ^ we present the ^-dependence of the second 
derivative d'^eo{0) /dO"^ on finite systems iV = 8, 10, 24, 26, 
as it follows from ( |4.2| ) . There is no signature in the finit e 
system results, which give a hint to the singularity (4.10). 
A comparison of (4.2) with a numerical differentiation of 
the ground state energy yields excellent agreement. 



FIG. 4. The second derivative of t he g round state energy 
eo{6) versus 9 as it follows from Eq. (4.2). 



V. LOW ENERGY EXCITATIONS IN THE SPIN 
SECTOR OF THE ID ANTIFERROMAGNETIC 
HEISENBERG MODEL 

The singular behavior of the first derivative of the stag- 
gered magnetization and of the second derivative of the 
ground stat e ene rgy per site can be traced back to the 
divergence ( 3.1C| ) of the quantity: 



/•OC 7 

E(6l = l,7r,7r,7V) = / — 5(w,7r,Ar), (5.1) 

which again can be viewed as the susceptibility of the 
dynamic structure factor 



S{cv,p, N) = Y, Mn{l,p)^S{u; - lu^). 



(5.2) 



The operators S(— p) • S(p), which enter in the definition 
( |3.4|) , only allow for transitions with = and Ap = 0. 
The gap uji{N) in the sector with quantum numbers of 
the ground state vanishes in the thermodynamic limit 



Nuji 25.15(5), 



(5.3) 



(4.9) 

as is shown in Fig. |M 



The divergence (3.10) of the susceptibility (5.1) de- 
mands for an infrared singularity in the dynamic struc- 
ture factor at p = TT 



(5.4) 



which is visible in th e nu merical data for the scaled tran- 
sition probabilities ( |3.4| ) of Fig. |[ in spite of the rather 
large finite-size effe cts. 

Insertion of (5.4) into (5.1) leads to ( 3.10| ). 

In Fig. 7 we present the p-dependence of the exci- 
tation energies Wn = En — Eq for = 24 together 
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FIG. 5. The excitation gap ujijN) for N = 8, 10, ... , 28. 
The solid line represents the fit (5.3). 



FIG. 7. The excitation spectrum for 9 = 1 and A'' = 24. 
The size of the dots is proportional to the relative spectral 
weights of associated transition matrix elements M„{l,p)^. 
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FIG. 6. The scaled transition probabilities (3.4) for finite 
systems (A'^ = 10, 28) vs. the excitation energy. 



with the associated relative spectral weights in percent- 
age terms. Highest spectral weights occur at excitation 
energies uJmax (p) ■ This curve is rather well approximated 

by 



i^maxip) = TTSinp. 



(5.5) 



The dashed curve is the dispersion relation for the lowest 
excitations in the St = 0-sector in the thermodynamic 
limit: 



■ smp, 



(5.6) 



as it was calculated in Ref. |T^. On our finite system we 
find nonvanishing excitation also below ujl(jj), they will 
not survive in the thermodynamic limit. 

Finally let us discuss the connected correlators: 

C{p,N) = {^o{l)\S{pfS{pf\^o{l)) 



^o{i)\s{pf\M^)r, 

dujS{cj,p,N), 



(5.7) 



whi ch fo llow by integrating up the d yna mic structure fac- 
tor (5.4). The infrared singularity (5.4) induces a diver- 
gence in C{tt, N) for N oo: 



(5.8) 



This behavior is clearly seen in Figjs. The slope of the 
numerical data in the inset o f Fig . ||is 0.52 5(8), which 
deviates from the prediction ( 3.11 ) and (|5.> ' 



A 



= 0.682(5). 



(5.9) 



by about -23%. 

The p-dependence of the static structure factor (5.7) 
is given in Fig. ^. One observes a singularity at p = n. 
Finite-size effects are small, away from the singularity. 
Near the singularity they can be desc ribed with a finite- 
size scaling ansatz in the fashion ( |3.14 ) we used for 
I]i(6' = l,7r,pi,7V) in Sec. III. 

In summary: The iV-dependence of the susceptibil- 
ity S(0 — l,7r,7r,iV) and of the connected correlators 
C(7r, N) as well as the p-dependence of 'S{9 = l,p,p, N) 
and C {p, N ) support the existence of the infrared singu- 
larity (5.4) in the dynamic structure factor (|5.2|). 



VI. THE ID ANTIFERROMAGNETIC 
HEISENBERG MODEL WITH NEXT TO 
NEAREST-NEIGHBOR COUPLINGS 

The approach developed in Sees. II - IV is applicable 
to all Hamiltonians H(i?) which depend on a parameter 
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VII. SUMMARY AND CONCLUSIONS 
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FIG. 8. The connected correlators (^) for iV = 8, . . . , 28 
are shown as sohd symbols (•), whereas the prediction of a 
finite-size scaling ansatz is marked by open symbols (o). The 
inset represents the singul ar b ehavior of C{n,N) together 
with the corresponding fit (|5.< 



■d and reduce to the ID Heisenberg Hamiltonian for the 
specific value = 1. As a further example we consider 
the Hamiltonian 



H(t?) = Hi + (i9 - 1)H2, 
with next to nearest-neighbor couplings 



N 



^ S„ • S 



n+2- 



(6.1) 



(6.2) 



71=1 



These couplings strengthen the antiferromagnetism for 
■d < 1, but frustrate the system for > 1. Therefore, 
the question arises, what happens here with the stag- 
gered magnetiz atio n for i9 < 1 . Repeating the arguments 
which lead to (3/7), yields for the first derivative of the 
staggered magnetization squared: 



= -2 



i5 = l 



it{1-2/N) 



dp 



cos(2p) -I- cosp] 



xE(i9 = l,p,TT,N). 



(6.3) 



Again I](z? = tt, iV) is given by ( |3.8[ ) in terms of the 
excitation energies LUmi'& = 1) = Ejn{'& — 1) — -Eo(i^ = 1) 
and transition probabilities Af,„(l,p) of the ID nearest- 
neighbor model with Hamiltonian H('!9 = JJ- However, 
insertion of the pole term contribution ( [3.12 ) leads here 
to a finite contribution for dm^ {d)'^ / d'd\i}^i in the ther- 
modynamic limit. It is worthwhile to note how the r ange 
of the spin-spin couphngs - over 2 lattice spacings in (6.1) 
=F 1 lattice spacings in (|2 



and over k = ^/N ^1 lattice spacings in ( |2.2| ) - enters 
via the Fo urier f acto rs cos(2p) and cos(fcp) in the for- 
mulas ( |6.3D and (3/7) for dm^dY / d'd\^^i, respectively. 
These Fourier factors generate th e di fferent beh avior of 
drn< {dY / d^U=i in the two cases (O) and (pj). 



It is convenient to study the transition from one to 
two dimensions by means of an interpolating Hamil- 
tonian H(6') [cf.( |2.2[ )] with a parameter 9, which con- 
trols the anisotropy between the horizontal and vertical 
nearest- neighbor couplings. The quantity of interest is 
the staggered magnetization m^{6), which on one hand 
is supposed to be nonzero for the anisotropic 2D model 
{Q <9 <1) and which on the other hand is known to be 
zero in the ID limit 0=1. In this paper we have studied 
the first derivative dm^{9Y/d9. It can be expressed in 
terms of excitation energies w„ = £"„ — Eq and transition 
matrix elements (^m|S(— p)S(p)|\E'o) associated with the 
operator S{—p) ■ S{p). 

In other words, the 0-evolution of the staggered mag- 
netization (|3.7|) is determined by the dynamics in the 
sector with total spin St = and momentum trans- 
fer Ap = 0. Energy conservation sum rules (3.(;) yield 



important constraints on the transition matrix elements 
(*,„(6l)|S(-p)S(p)|*o(6')). They turned out to be very 
helpful in our numerical analysis of the low-energy exci- 
tations at 6* = 1: 



1. The dynami c st ructure factor ( |5.2D has an infrared 
singularity (5.4) at p = tt. For < p < it, there is 
a dispersion relation uj = Wi(p) ( [5.6[ ) for the low- 
est excitations (with AS* = and Ap = 0) in the 
thermodynamic limit. On finite systems, we find 
excitations below lul (p) ■ 

2. The existence of the infrared singularity is sup- 
ported by the divergences, which can be clearly seen 
in the N dependence of the connected correlators 
C{tt,N ) (U) and the susceptibility S(l, tt, tt, A^) 
( 3.10 ) and in the p-dependence of C{p,N) and 
S](l,p,p, N) for p ^ TT. 

3. Owing to the energy conservation sum rule ( [3.S| ), 
t he d ivergence in the susceptibility induces a pole 
( 3.12 ) in S(l,p, tt) . Th is pole term fixes the leading 
large- A behavior (3.15) of the first derivative in the 
staggered magnetization squared at 9 = I. 

We have considered in this paper the transition from 
one to two dimensions in the spin 1/2-antiferromagnetic 
Heis enberg model o nly . However, the general formulas 
( |3.3[ ) and sum rule ([3.6|), which describe the 6'-evolution, 
hold as well in the spin-1 case, where the IXli-Lnodel 
(9 = 1) develops the celebrated Haldane gap.t^^3 The 
existence of this gap already t ell us, that there is no in- 
frared singularity of the type (gj). Therefore, we do not 
expect a divergence (3.15) in dijv {9Y / d9\e=i in the spin 
1 case 



APPENDIX: RECURSION METHOD 

The recursion method is designed to approximately de- 
termine the excitation energies and matrix elements of 
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transition operators, the dynamics of which are to be 
studied. Neglecting indices Uke momentum p and chain 
length N let us consider 

fAico) = J2 l("l^|0)l''5[c. - - Eo)], (Al) 



with |0), \n) being the ground state and the excited states 
of the system defined by the Hamiltonian H. Eq and _E„ 
are the corresponding energies, (e.g. A might be chosen 
a.s A — Sj{p,N), j — 1,2,3). The Laplace transform 
of //I (a;) reads 



/^(r)=Do(T)(/o|/o) 



(A2) 



{0\A+t 



^|0) = (/o|/(r)), (A3) 



with l/o) = A\0) and |/(t)) fulfilhng 

|-|/(r)) = -H|/(r)), (H = H - Eo). (A4) 

Now, a Gram-Schmidt construction is used (see also Ref. 
^ ) to form an orthogonal set of states {l/fc)}, namely 



|/,+i)=H|/,)-afc|A.)-6^|/;c-i) 



(/fc|H|/fe) 

gfc = , , I ^ , , fe = 0, 1,2, ..,L - 1, 



&o = 0, 

bl = , y^^/f' , , A: = 1,2,..,L- 1, 



(/fc-il/fc-i) 
and |/(t)) is expressed as 

L-l 



|/(r)) = 5]D,(r)|/fc). 



(A5) 

(A6a) 
(A6b) 
(A6c) 

(A7) 



fe=0 



L denotes the dimension of the Hilbert space, i.e. the 
number of states \n) leading to no nzer o matrix elements 
(7i|A|0). The differential equation ( |A^) leads to 

^ Dfe(T)|/fe) - - ^Dfe(r) (lA+i) + ak\h) + bl\h-i)) 

k k 

(A8) 

or in matrix notation D = — MD with M being a tridi- 
agonal L x i-matrix with eigenvalues a;„ and eigenvectors 
en- The eigensolutions Vkij) of the set of linear differen- 
tial equations read Wfc(r) = w°e~"'''^ leading to 



Dfc(r) = ^{en)kv\ 



(A9) 



Here, the are given by the initial condition D(r = 0) = 
(1,0,0,...,0)^ yielding for /^(r) 



^'°(e„)o(/o|/o), (AlO) 



i.e. |(n|^|0)|2 = u°(e„)o(/ol/o)- It should be noted that 
the eigenvalues of M are real, however, its non non- 
symmetric structure does not require the e„ to be or- 
thogonal. The approximation in our scheme sets in in 
the fact that the numerical treatment does not allow for 
generating the whole se t o f L orthogonal states. The 
number of iterations in (A5) is reduced to L, the matrix 
M truncated to L x L and the evaluation given above 
is performed on this reduced set of states. Details and 
checks for this method are given in Refs. pT| , p2[ 

In the present paper the recursion method has been 
extended [see e.g. Eq. ( |3.3| )] to treat matrix elements of 
type 



(5|/(r)) = (.9|e-"|/o), 



(All) 



with l/o) = A|0), \g) = B|0), A and B acting on the 
same space. Using the series representation of |/(r)) we 
obtain 



(0|B+|n)(n|A|0)=^u°(e„)fe(5|/. 



(A12) 



which contains the transition probabilities |(n|A|0)p for 
the choice B = A. 
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